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This pedagogical comment highlights three misconceptions concerning the usefulness of the con¬ 
cept of negative temperature; being derived from the usnal, often termed Boltzmann, definition 
of entropy. First, both the Boltzmann and Gibbs entropies mnst obey the same thermodynamic 
consistency relation. Second, the Boltzmann entropy does obey the second law of thermodynamics. 

Third, there exists an integrating factor of the heat differential with both definitions of entropy. 


I. INTRODUCTION 

The concept of negative temperature first arose, for 
spin systems, in the 19509^*^ and has since been useful 
for understanding a variety of systems with spontaneous 
magnetic order—. The essential feature of these systems is 
an energy spectrum that is bounded above, as well as be¬ 
low. This means that the canonical partition function can 
undergo analytic continuation to negative temperatures. 
The total energy of the system at negative temperature 
is higher than that at infinite temperature. 

More recently, negative temperatures have been ob¬ 
served relating to motional degrees of freedom^. An 
ensemble of particles with an attractive interaction, or 
negative pressure, was prepared confined to a band with 
upper and lower bounds on kinetic energy. Since the sys¬ 
tem was observed to be stable, not imploding, we can 
conclude that its temperature must be negative. 

Since the idea of temperatures being negative is unin¬ 
tuitive, there has arisen a controversy as to whether it 
may be advantageous to use an alternative definition of 
entr^y that guarantees only positive temperatures. In 
ref. yit has been suggested that we should use the Gibbs, 
or volume, entropy since it increases monotonically with 
energy. 

To elaborate we consider an isolated system with en¬ 
ergy E in the microcanonical ensemble^. We denote as 
T{E) the number of microstates that can be accessed by 
the system. It is important to recognise that r(i?) is a 
distribution. Especially if the energy spectrum is contin¬ 
uous then it gives the number of microstates in a distri¬ 
bution bin of some width A. The cumulative distribution 
Ul{E) represents the sum of T{E) from the ground state 
up to energy E; we assume that the energy spectrum is 
bounded below. 

The two definitions of entropy are the Boltzmann, or 
surface, entropy 

SB = kB\nT{E) ( 1 ) 

and the Gibbs, or volume, entropy 


tially. In these circumstances, both entropy definitions 
are expected to provide the same results^. 

The temperature T of the isolated system is given by 


Since the cumulative distribution n{E) is uniformly in¬ 
creasing, the Gibbs temperature is necessarily positive; 
Tg > 0. In contrast, for systems with an energy spec¬ 
trum bound above, it is the case that T{E) is a decreasing 
function for high energy and this means that the Boltz¬ 
mann temperature Tb can be negative. 

The main proposal of ref. is that it is always appro¬ 
priate to use the Gibbs entropy and reject the Boltzmann 
entropy. This has created a controversy with support for 
the proposali^— as well as criticism in support of keeping 
the Boltzmann entropy and negative temperatures^i^— . 
In ref. [Ill, for instance, it is shown that using the Gibbs 
entropy can mean heat flowing from cold to hot as well 
as the temperature Tq not being intensive. 

The aim of this short note is not to review the con¬ 
troversy in detail but just to point out three simple ob¬ 
servations in support of negative temperatures and thus 
provide a clean bill of health for the Boltzmann entropy. 
First, it is not true that thermodynamic consistency is 
ever violated if we use Sb- Second, use oi Sb does not 
violate the second law of thermodynamics. Third, we can 
find an integrating factor for the heat differential using 
either definition of entropy. 


II. THERMODYNAMIC CONSISTENCY 

Suppose we have three variables S, E and h and that 
some function of them obeys f{S, E, h) = 0. Then 

df = fsdS + fsdE + fhdh = 0 (4) 

where the coefficients are partial derivatives. First, fixing 
h, we can prove that 


SG = kBlnn{E) ( 2 ) 

where kB is Boltzmann’s constant. The usual situation 
in statistical mechanics is that the energy spectrum is not 
bounded above and the value of r(i?) increases exponen¬ 






( 5 ) 


Then, fixing each variable in turn, we derive three equa¬ 
tions from which we can eliminate the derivatives of / 


2 


and arrive at 


to deduce that 






( 6 ) 


These equations are very well known from 
thermodynamics^. With T defined by 


1 

T 



(7) 


we can see that, from dH]), 


T \^)s 


-1 


( 8 ) 


E <{H) (15) 

It is important that all that has physical significance is 
independent of the bin width A. Although we choose A 
to be very small it should not depend on the number of 
particles N. For example A = e~^ would not provide 
us with a distribution. We must insist that the total 
number of distribution bins n is very small in relation to 
the number of particles; that is n N. Now, since the 
energy E is an extensive quantity, we have that E ^ N 
and can treat A <§; A in the thermodynamic limit^. Thus 

E = {H) (16) 


and then, using a variant of ([5]), 


T 




(9) 


So far this is just mathematics; S{E,h) can be any 
single-valued section of a differentiable function without 
singular points. We continue by considering an isolated 
system with entropy S', fixed energy E and subject to an 
external influence h. Thermodynamic consistency says 
that^ 



where the El is the Hamiltonian and the average is eval¬ 
uated in the microcanonical ensemble. The influence h 
can of course be interpreted as a vector if there are many 
components. 

First, let us consider a classical system with Hamilto¬ 
nian iJ (p, q). As usual, p and q are vectors for the canoni¬ 
cal momenta and coordinates. In three space dimensions 
with N particles we have vectors in a 3A-dimensional 
space. The energy spectrum is not discrete and we need 
to identify F(i?) as the number of accessible microstates 
in a distribution bin of some width A. We can write^ 

F(i?) = j dpdq (11) 

J E<H<E+A 


with sufficiently large N. 

Returning to consider Eq. (nni we have that 

We can imagine a second microcanonical ensemble with 
the same entropy S but subject to an influence h + dh. 
The energy of this second ensemble must he E + dE with 

The way to fix the entropy in Eq. (II3 is to simply inte¬ 
grate over the microstates accessible to the first unper¬ 
turbed ensemble. Then the derivative only acts on the 
integrand and 

If dH ^ ^ 

\dhjg h{E) joh 

and this is precisely what we desire. Of course, fixing 
the set of accessible microstates also fixes T{E) and the 
entropy Sb = fcslnF(A). Furthermore we have also 
fixed the cumulative distribution El(E) and the entropy 
Sg = ks In 0(A). It does not matter which definition of 
entropy we use to prove thermodynamic consistency. 

It is worthwhile to emphasize here that we are just 
using the obvious statistical expression 


In order to make this dimensionless it is usual to divide 
by Planck’s constant 3N times but this detail will not 
prove important. 

The ensemble average of the Hamiltonian is 


(H) 


^Je<h<e+a 

^ \^) JE<H<E+A 


( 12 ) 

(13) 


Then we can use the bound 


1 

T{E) 


IE<H<E+A 


{Ed — E)dpdq 


< A 


(14) 



where the averages over the microstates must have ex¬ 
actly the same meaning. The mapping from a set of 
microstates to some thermodynamic entropy function is 
essentially free, although perhaps unphysical in general. 

For systems with a discrete energy spectrum, it is most 
convenient to put one energy level into one distribution 
bin, although we might in principle put more than one. In 
any case we just need to sum over the levels in the bin. 
For spin systems we also need to replace the integrals 
over p and q with discrete sums. However, regardless 
of these details, we fix the entropy by a trace over the 









3 


microstates accessible to the unperturbed system. We 
have thermodynamic consistency for all systems using 
either Sb or Sq- 

This can be illustrated by the example of a gas of Ising 
spins. The Hamiltonian is 


N 

H = -h'^s, ( 21 ) 

i=l 


Here h is an external field with units of energy. The 
Ising spins take values Si = ±1. This is basically the 
two-level system commonly introduced to illustrate the 
idea of negative temperature^i^ d^d^d^i^^d’^ . The energy 
spectrum is discrete; with up spins with Si = -1-1 and 
Ni down spins 


E={H) = -h{N^ - N^) (22) 


Since the spectrum is discrete we put one energy level 
into one distribution bin and perform discrete sums over 
the spins. With up spins the number of accessible 
microstates with energy E is 


r{E) = 


Nl 


(23) 


where and N = N^+N]^. The cumulative 

distribution is 


Nx 


m) = E 


TV! 


n=0 


n\{N — n)! 


(24) 


As we expect, fixing the value of fixes r(i?), Sb, H(i?) 
and Sg- Our thermodynamic consistency equation reads 


T 





where we can use either definition of entropy. 
Our general conclusion is that 


Tb 




(26) 


III. SECOND LAW OF THERMODYNAMICS 

We aim to prove here that the Boltzmann entropy Sb 
is always consistent with the second law of thermody¬ 
namics. In ref. 0 it is suggested that this may not be the 
case under certain conditions. Following ref. 0 we state 
the second law, according to Planck, as 

S{Ex + E 2 ) > Si{Ex) + S2{E2) (28) 

The terms on the right are the entropies of two mi- 
crocanonical ensembles with respective energies Ei and 
E 2 ■ We imagine coupling these systems and waiting long 
enough for a new thermodynamic equilibrium. Then the 
term on the left is the entropy of the composite system 
which has energy Ei+ E 2 . The equality occurs if the un¬ 
coupled systems have the same temperature. Otherwise 
the entropy increases. 

The uncoupled systems have distributions F 1 {Ei ) and 
T 2 {E 2 ) for their numbers of accessible microstates. Thus, 
before coupling, the number of microstates accessible to 
the composite system is Ti{Ei)T 2 {E 2 ). After coupling, 
the distribution for the composite system will be^ 

n 

r{Ei +E2)=Y, Ti{E,)r 2 {Ei +E 2 -E,) (29) 

i=l 

with a sum over n ^ N distribution bins. The uncoupled 
systems have numbers of particles proportional to N and 
their distributions use the same numbers of bins. On 
the left we should best use a bin width 2 A double that 
on the right. With a discrete spectrum we can take a 
simple sum over the energy levels. Otherwise, as long as 
we maintain n <C TV, we can take n to be very large and 
write 


pEi -\-E2 —Eo2 

riEi+E 2 )= Ti{E)r 2 {Ei+E 2 -E)dE 

J Eq\ 

(30) 

where Aqi and A 02 are the ground state energies of the 
uncoupled systems. 

In all cases we have that 


in contradiction to a statement given in ref. 0. A proof 
that the Boltzmann entropy Sb does obey thermody¬ 
namic consistency has also been given in ref. [T^. The 
technique included the saddle-point method which is cer¬ 
tainly valid where N is sufficiently large. For the gas of 
Ising spins the energy, or magnetic enthalpy, is simply 

E = —Nhta,nhf3h (27) 


r{Ei + E2)>ri{Ei)r2{E2) (3i) 

since all the terms in the sum are nonnegative. If we 
now define entropy according to S = fcs In F then we 
arrive at the inequality (l2^ . This is a proof that using 
the Boltzmann entropy never violates the second law of 
thermodynamics. In fact, the sum (1291) is dominated by 
one particular term^ 


where (3 = and the temperature T is interpreted 
as the Boltzmann temperature Tb- If more than half of 
the spins are down then the energy is positive and the 
temperature must be negative. This result can be de¬ 
rived simply from either the canonical or microcanonical 
ensemble. 


r{Ei+E2)-r,{Ei)T2{E2) (32) 

where Ei and E 2 are the energies corresponding to ther¬ 
mal equilibrium. This leads to 

SiEi+E2) = Si(Ei) + S2(E2) (33) 
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in the thermodynamic limit; in agreement with the pos¬ 
tulate of the additivity of entrop} d'^d^ which is essential 
with regard to prediction of the equilibrium energies. 

The density of states oj{E) is defined^ by 

r(£;) = ui{E)A (34) 

Provided the distribution bin width A is sufficiently 
small, we can write uj{E) = ^ . In ref. 0 it is written 

that r = ecu where e is described as a parameter that 
ensures that P is dimensionless. In fact, e must be inter¬ 
preted as a distribution bin width; not freely. In terms 
of the density of states, inequality dSD reads 

u}{Ei + E 2 ) > cui(i?i)cu2(i?2)—^— (35) 

We can define^ entropy via S = fcs Incu and, with say 
Ai = A 2 = iA, 

S{Ei+E 2) > Si{Ei) + S2iE2) + kB\n^A (36) 

The third term must be negligible in the thermodynamic 
limit since entropy is extensive, S N. As remarked 
above, choosing say A = e~^ is not possible since we 
would not have a distribution for r(i?). The value of 
A must not depend on N and both S = Acs In P and 
S = ksEiui obey the second law. 

Much discussion as to whether the Gibbs en¬ 
tropy ks In Q obeys the second law can be found 
elsewhereiiiid^dii^— . In ref. 0 it is proven that 

n{Ei + E 2 ) > ni{Ei)n2{E2) (37) 

and we can conclude that the Gibbs entropy Sc does sat¬ 
isfy the second law as stated in Eq. (E51) . Nevertheless, 
the crucial issue concerns the additivity of entropy as ex¬ 
pressed in Eq. In the event that the Boltzmann 

entropy in increasing with energy it is well knowni^^ that 
we can prove that Sb and Sq provide the same ther¬ 
modynamic predictions. It is in the case that Sb is not 
increasing that difficulties aris o^'^d^d'^ . 

Following the construction in ref. 0 we can write ex¬ 
actly 

n(Ei + E2) = ni(Ei)n2(E2) 

PE2 — E02 pE 2 — E 

+ / Ti{E + Ei) / T2{E')dE'dE 

Jo J Eq 2 

f'E\ — Eq\ pEi — E 

+ / T2{E + E2) / Ti{E')dE'dE 

Jo J Eqi 

(38) 

The Gibbs entropy will only give a correct thermody¬ 
namic prediction if the second and third lines do not 
grow exponentially with system size. Looking at the sec¬ 
ond line above, we can image that T 2 {E') is not mono¬ 
tonic and possesses a peak at some value E' = e' . If 
E < E 2 then the line will indeed provide a contribu¬ 
tion that grows exponentially with size. This destroys 
the utility of the Gibbs entropy. 


To illustrate we can use the gas of Ising spins. It is con¬ 
venient here to parameterise the energy with the number 
IV 4 , of down spins; the ground state has = 0. We cou¬ 
ple two systems with Ni and N 2 spins where N = Ai-|-iV 2 
and Ni ^ N 2 ^ N. The total number of down spins 
= N^i -|- is conserved. In thermal equilibrium we 
have 


N II _ N ^2 _ 

Ni ~ N2 ~ N 

Translating the second line of (1381) . we have 


(39) 


N ^2 Ni_ 2—'^1 

^ ri(TOi -I-iV4,i) ^ T2{m2) 


m\ — l 


m2—0 


(40) 


where 


r 2 (TO 2 ) 


N 2 I 

77121(^2 — m 2 )! 


^N2f(x) 


(41) 


with X = ^ and f{x) = —x\nx — (1 — a;)ln(l — x). It 
is a simple matter to show that f{x) has a maximum at 
cc = d, or m 2 = \Al 2 - This will contribute something 
large to the sum (gOl) if N ^2 > \N 2 or, equivalently, 
Ni > dfV which is precisely where the temperature is 
negative and the Gibbs entropy fails to provide any ther¬ 
modynamic prediction of the equilibrium state. A nu¬ 
merical evaluation of this has been given in ref. [13 where 
it is also clear that the Gibbs temperature Tc is not in¬ 
tensive, depending on system size contrary to common 
sense. 


IV. INTEGRATING FACTOR FOR THE HEAT 
DIFFERENTIAL 


In thermodynamics the heat differential dQ is not 
exact^. The integrating factor y, where T is temper¬ 
ature, allows us to use an exact differential 


dS='^ (42) 

where S is entropy. In ref. 0 it is suggested that this 
might not be valid in every case unless we use the Gibbs 
entropy Sc and the Gibbs temperature Tq. 

The work done, on a system, by an external influence 
h is 

where we follow the notation above. Thus 

dQ = dE-dW = dE + T dh (44) 

using Eq. (US. We know that we can use either the 
Boltzmann or Gibbs definitions of entropy, and temper¬ 
ature, here. 
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Suppose we use S = the Gibbs definition. 

Then, since T = (§f )^, 

dQ = dE+-^dh (45) 

oj ah 

where a; = ^ is the density of states at fixed h and the 
derivative of ^{E, h) is at fixed energy. We can see that 
u! is an integrating factor since 

BQ BQ 

ojdQ = -T—dE + -T^dh = did (46) 

BE Bh 

Continuing, ref. 0 correctly argues that integrating fac¬ 
tors for dQ must be of the form where g is differ¬ 

entiable. The law of ideal gases then identifies a unique 
solution g{il) = /cb Infl. 

The point now is that we can return to Eq. (PI) and 
use the definition S' = As In T for the entropy and simply 
replace il with T in all expressions. We can do this since 
both definitions obey the same thermodynamic consis¬ 
tency relations; we refer to Eq. (1261) . In summary, both 
entropy definitions provide an integrating factor 

As a matter of fact we can replace il with any S as 
long as, in the case of the ideal gas, 

EP (47) 

with p = ^ — c where c is constant, or unimportant in 
the thermodynamic limit. We have c = 0 with il and 
c = 1 with r. 


V. DISCUSSION 

We have proven three important properties of the 
Boltzmann entropy that have been challenged in the lit¬ 
erature. It always obeys thermodynamic consistency, the 
second law of thermodynamics and provides an integrat¬ 
ing factor for the heat differential. 

One important theme is that the number of accessible 
microstates r(£’) should be regarded as a distribution. 
It does not necessarily have to be differential or even 
continuous although we usually assume that it is. 

A possible further comment is that we have the mi- 
crocanonical ensembles subject to an external influence 
when they are supposed to be isolated. Also, for the gas 
of Ising spins, we have no interactions between the parti¬ 
cles and may wonder how two systems can mix or couple 
and find a new thermal equilibrium. The usual answer 
to this is that the particles have weak interactions. 

These points might be addressed by looking at the ex¬ 
ample of an Ising model in one dimension with Hamilto¬ 
nian 

N-l 

H = SjSi+i (48) 

2=1 


which has also proven usefnl in ref. H^ . The exchange in¬ 
teraction J is positive, as large as we like, and there is no 
external influence. To solve this problem we can imag¬ 
ine a bond with parallel spins as represented by another 
Ising variable cr = -|-1. A bond with spins not parallel 
has tr = — 1. We take N to be very large and do not 
worry about boundary effects. 

The degeneracy of the energy level E = —J{N^ — N_) 
is given by 

The factor of 2 is for global inversion of the spins. 
is the number of bonds with parallel spins and N- = 
N — N+. In the same way as for the gas of Ising spins we 
can solve this model by just writing J instead of h. We 
have that the total energy is 


E = -NJtanh^J (50) 


If more than one half of the bonds have spins not parallel 
then this energy is positive and the temperature is nega¬ 
tive, provided we use the usual Boltzmann definition. 

We can imagine some divine intervention that arranges 
a majority of bonds with spins not parallel and then iso¬ 
lates the system very quickly. The spins would then come 
to thermal equilibrium, through their interactions, at a 
negative temperature. 

A similar illustration can easily be considered in higher 
dimensions. For instance, a square lattice Ising model, 
with all exchange interactions J positive, can be isolated 
at negative temperature if a majority of bonds have spins 
not parallel. A perfect antiferromagnetic arrangement 
would indicate the extreme case of temperature negative 
zero. 

An important utility of the concept of negative temper¬ 
ature is that we can use the canonical ensemble. Where 
the energy spectrum is bounded above, the partition 
function^ can be defined at negative temperature by sim¬ 
ple analytical continuation. Since the microcanonical en¬ 
semble is generally much less mathematically tractable, 
this is a powerful point. 

The Gibbs temperature is guaranteed to be positive 
and this is intuitively attractive. Nevertheless, we cannot 
use the canonical ensemble and there are serious issues 
about thermodynamic prediction in the case that r(i?) 
is not monotonically increasing . 
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